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We investigate the structural aspects of genuine multipartite entanglement in Gaussian states of continuous 
variable systems. Generalizing the results of [Adesso & Illuminati, Phys. Rev. Lett. 99, 150501 (2007)], 
we analyze whether the entanglement shared by blocks of modes distributes according to a strong monogamy 
law. This property, once established, allows to quantify the genuine A^-partite entanglement not encoded into 
2, . . . , iC, . . . , (A^ — 1) -partite quantum correlations. Strong monogamy is numerically verified, and the explicit 
expression of the measure of residual genuine multipartite entanglement is analytically derived, by a recursive 
formula, for a subclass of Gaussian states. These are fully symmetric (permutation-invariant) states that are 
multi-partitioned into blocks, each consisting of an arbitrarily assigned number of modes. We compute the 
genuine multipartite entanglement shared by the blocks of modes and investigate its scaling properties with the 
number and size of the blocks, the total number of modes, the global mixedness of the state, and the squeezed re- 
sources needed for state engineering. To achieve the exact computation of the block entanglement we introduce 
and prove a general result of symplectic analysis: Correlations among K blocks in A/^-mode multi- symmetric 
and multipartite Gaussian states, that are locally invariant under permutation of modes within each block, can 
be transformed by a local (with respect to the partition) unitary operation into correlations shared by K single 
modes, one per block, in effective non-symmetric states where N — K modes are completely uncorrelated. Due 
to this theorem, the above results, such as the derivation of the explicit expression for the residual multipartite 
entanglement, its non-negativity, and its scaling properties, extend to the subclass of non-symmetric Gaussian 
states that are obtained by the unitary localization of the multipartite entanglement of symmetric states. These 
findings provide strong numerical evidence that the distributed Gaussian entanglement is strongly monogamous 
under and possibly beyond specific symmetry constraints, and that the residual continuous-variable tangle is 
a proper measure of genuine multipartite entanglement for permutation-invariant Gaussian states under any 
multi-partition of the modes. 

PACS numbers: 03.67.Mn, 03.65.Ud 



I. INTRODUCTION 

Although multipartite entanglement is ubiquitous in the 
physics of many-body quantum systems, it is very difficult 
both to characterize and to quantify it (ll 0]. While a rather 
satisfactory understanding has been achieved in the bipartite 
case 1 3, 4], there is a certain degree of consensus that there 
is no unique way to define multipartite entanglement, even in 
the simplest case of pure states L3|, |5J . The presence of mul- 
tipartite entanglement clearly depends on the partitioning that 
one imposes in order to group the individual subsystems into 
parties. Furthermore, given a fixed partition one can single out 
a hierarchy of different levels of multipartite entanglement, or 
alternatively /c- separability, which establishes a smooth con- 
nection between the two limiting cases of fully separable state, 
where the parties are all disentangled, and of fully insepara- 
ble states, where entanglement exists across any global bisec- 
tion, and the parties are supposed to share genuine multipar- 
tite entanglement. To complicate things further, not even the 
concept itself of full inseparability is uniquely defined, as in 
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all quantum systems (starting with the simplest paradigmatic 
one: a system of three qubits [6]) there exist several inequiv- 
alent forms and classes of fully inseparable states, with dif- 
ferent properties and applications. Therefore, different routes 
to the definition, characterization, detection, and quantifica- 
tion of multipartite entanglement are possible and have been 
in fact proposed during the last decade yj,|3|,l5|]. 

In this context, we would like to recall here three relevant 
approaches, that are of interest from different perspectives. 
One possible way, directly generalizing the concept of bipar- 
tite entanglement, is to quantify the multipartite entanglement 
present in a quantum state as the distance between the state 
itself and the set of fully separable, or (more finely-grained) 
/c-separable states, thus adopting an axiomatic approach and a 
geometric definition of entangled states in terms of their "dis- 
tinguishability" from separable states |7, 8]. Another possi- 
bility is to regard entanglement as a physical resource and ap- 
proach its definition operationally by relating different mea- 
sures of multipartite and global entanglement to the success 
or the performance of various multi party quantum informa- 
tion and communication protocols iSllSllll]- A further pos- 
sibility stems from the observation that bipartite entanglement 
distributes obeying a so-called "monogamy" law 1 12]. From 
the mathematical formulation of such a property, constraints 
are imposed on genuine multipartite entanglement that lead, 
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in special cases, to suitably defined measures for its quantifi- 
cation (H [Hill [H. 

The above, together with other definitions such as quantita- 
tive entanglement witnesses [17] and axiomatic measures like 
the multipartite squashed entanglement fT§], all come with 
their burden of mathematical difficulties, e. g. in the solu- 
tion of the involved optimization problems, and/or with their 
physical and operational limitations, e.g. in the dimension- 
ality of the systems to which they are applicable, or in the 
specific protocol they rely on. The various approaches that 
we have briefly reviewed have nevertheless enabled important 
progress in the quest for the characterization of multipartite 
entanglement. Most of these progresses have been achieved 
for systems of qubits, as the complexity in the description of 
quantum states scales exponentially with the Hilbert space di- 
mensionality of the constituents. However, the situation be- 
comes more tractable if one moves to algebras of observables 
with continuous variable (CV) spectra in infinite-dimensional 
Hilbert spaces [19, 20], and in particular if the analysis is re- 
stricted to multimode Gaussian states. In this case, one can 
derive directly some basic results in analogy with systems of 
qubits; moreover, some very interesting new features emerge 
1I21I1 . The obvious advantage, when dealing with Gaussian 
states (which include states of fundamental physical inter- 
est, such as vacuum, coherent, squeezed, and thermal states 
of harmonic systems) is that they are entirely characterized 
by the first and second moments of the canonical quadrature 
phase operators, which delimits the object of the investigation 
to a finite set of degrees of freedom. Due to this drastic di- 
mensional reduction and to their physical relevance, Gaussian 
states are still, at present, the theoretical and experimental pil- 
lars of quantum information with CVs 12211 . 

Concerning the multipartite entanglement of Gaussian 
states, the last two approaches out of the three mentioned 
above have been particularly successful. From an opera- 
tional perspective, important basic results were obtained by 
van Loock and Braunstein [23]. These authors have intro- 
duced schemes to produce fully symmetric (i.e. permutation- 
invariant) N-mode Gaussian states exhibiting genuine multi- 
partite entanglement and, accordingly, have devised a com- 
munication protocol, the "quantum teleportation network", 
for the distribution of quantum states exploiting such re- 
sources. Subsequently, we have elucidated that the optimal 
nonclassical fidelity characterizing the performance of such 
a protocol (which has been experimentally demonstrated for 
N = 3 [24]) is in one-to-one correspondence with the pres- 
ence of genuine multipartite entanglement in fully symmet- 
ric Gaussian resources [10]. This operational characterization 
of Gaussian entanglement in the framework of quantum tele- 
portation has been followed by a series of results establish- 
ing the ("weak'O monogamy of entanglement for all Gaus- 
sian states fT3, "26]. While in its traditional expression 
the monogamy of entanglement is known to constrain only 
the bipartite entanglement in different partitions of a multi- 
partite system (as originally proven for qubits (v3, '27]), in 
recent work we have defined a broader framework that pos- 
tulates the subsistence of a stronger monogamy inequality. 
This property of "strong monogamy" imposes the existence 



of a trade-off on bipartite and genuine multipartite entangle- 
ment simultaneously [28]. We have demonstrated the valid- 
ity of our construction for the entanglement shared by single 
modes in the important instance of fully symmetric (pure or 
mixed) Gaussian states of an arbitrary number N of modes, 
that constitute the set of entangled resources in the van Loock- 
Braunstein teleportation network. Remarkably, it is possible 
to establish that the information-theoretic measure of multi- 
partite entanglement emerging from such a strong monogamy 
decomposition — dubbed "residual contangle" — is a mono- 
tonic function of the optimal nonclassical fidelity of teleporta- 
tion for teleportation networks. This result integrates the op- 
erational and "distributional" approaches to the quantification 
of genuine multipartite entanglement when applied to multi- 
party Gaussian states with complete symmetry under permu- 
tation of the modes. In passing, we recall that other applica- 
tions of symmetric, multiparty entangled Gaussian states have 
been proposed and demonstrated, including secret state shar- 
ing [29] and Byzantine agreement |30]. More details on these 
and other applications of multipartite Gaussian states in quan- 
tum information can be found in Ref. lHO, [22I], as well as in 
two recent complementary reviews ll2lll (theoretical) and lI sTI] 
(experimental). 

In the present work we build on the ideas presented in Ref. 
[28] and consider the quantification of genuine multipartite 
entanglement of Gaussian states under less stringent sym- 
metry constraints and more general partitions of the modes, 
comprising the instance of an arbitrary number of modes 
per party. We first analyze in general the postulated strong 
monogamy decomposition: Specializing to fully symmetric 
Gaussian states, we derive an analytical formula for the resid- 
ual contangle. We then show that in such states the strong 
monogamy decomposition holds for the entanglement shared 
by arbitrary- sized blocks of modes (in other words, that the 
corresponding residual contangle is nonnegative), providing 
strong numerical evidence which extend the domain of valid- 
ity of the analytical proof of Ref. i28ll (valid for entanglement 
shared among single modes). This result yields a computable 
tool that can be readily applied to quantify genuine multipar- 
tite entanglement of symmetric Gaussian states. To this aim, 
we include some handy numerical and analytical-algebraic 
codes that can be readily computer-implemented. We then 
investigate in detail how the multipartite entanglement scales 
with the total number of modes, with the number of parties in 
which the global system is partitioned, with the size of the in- 
dividual blocks of modes pertaining to each party, and, finally, 
with the resource parameters such as the degree of squeezing 
applied on each mode. This analysis, corroborated by numer- 
ous explicit examples, allows to gain a finer understanding of 
multipartite Gaussian entanglement and to identify its charac- 
teristic traits, as well as to single out the optimal conditions in 
which multipartite entangled Gaussian states can be employed 
as practical resources for quantum communication tasks. 

Here and in the following we focus on Gaussian states of 
N modes, partitioned into K blocks, dubbed "molecules" 
for ease of visualization [28], under the constraint of multi- 
symmetry. Imposing the latter property means that we con- 
sider states invariant under permutations of modes within each 
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molecule. Clearly, fully symmetric states are those states that 
are multi- symmetric with respect to any multi-partition. Gen- 
eralizing a result proven for bipartite bi-symmetric Gaussian 
entanglement by Serafini, Adesso, and Illuminati [32], we 
prove that the "molecular" genuine i^T-partite Gaussian entan- 
glement shared by K blocks is unitarily localizable. Specifi- 
cally, this means that by performing only local unitary opera- 
tions independently in the subspace of each single molecule, 
multi- symmetric A^-mode Gaussian states are reducible to an 
equivalent form in which only K modes {K < N), one per 
molecule, are correlated — we call these target modes "nu- 
clei" — while the remaining {N — K) modes become com- 
pletely uncorrected. Therefore, the molecular entanglement 
in the original multi- symmetric A/'-mode Gaussian state across 
a given AT-block partition is equal to the nuclear entanglement 
in the localized Gaussian state of K modes. Interestingly, 
such a localized nuclear state does not exhibit, in general, 
any specific symmetry under permutations of the modes. This 
fact allows to test explicitly the strong monogamy decomposi- 
tion i28ll on a restricted class of random pure and mixed non- 
symmetric multimode Gaussian states as well, obtained by 
unitarily localizing entanglement in fully symmetric Gaussian 
states of a higher number of modes with a random structure of 
partitions. After an extensive performance of such a test, no 
violation of the strong monogamy constraint has been found. 
This result provides a first evidence that CV Gaussian entan- 
glement can be strongly monogamous even beyond specific 
symmetry requirements, and that the integrated operational- 
distributional approach to the quantification of CV entangle- 
ment might possibly lead to a proper measure of genuine mul- 
tipartite entanglement in general, fully inseparable multimode 
Gaussian states. On the other hand, more extended tests in 
more general instances are needed to fully validate such a 
strong conjecture. 

The plan of the paper is as follows. In Sec. [Ill we re- 
view some basic notions on Gaussian states of CV systems, 
some fundamental results on the characterization and quan- 
tification of Gaussian entanglement, and some relevant facts 
on distributed entanglement and monogamy constraints. In 
SecJIIIwe briefly discuss the route introduced in Ref. |28] to 
the definition of a measure of genuine multipartite entangle- 
ment in CV systems by postulating the existence of a stronger 
monogamy inequality that constrains both bipartite and mul- 
tipartite entanglement. In SecjlVlwe prove that correlations 
among blocks of modes in multipartite Gaussian states which 
are multi- symmetric with respect to local mode permutations 
can be completely and reversibly localized into effective cor- 
relations shared by single modes under a local unitary opera- 
tion. In SecOwe apply this result and discuss the quantifica- 
tion of genuine multipartite entanglement distributed among 
molecules of arbitrary size in globally pure or mixed, fully 
symmetric Gaussian states. We derive an analytical expres- 
sion for the entanglement and study its scaling properties with 
the relevant physical parameters, including purity, squeezing, 
number of modes, number and size of molecules, and the spe- 
cific multi-partitions imposed on a multimode CV system. In 
Sec.|Vllwe exploit the property of unitary localization to re- 
interpret our findings in terms of the genuine multipartite en- 



tanglement shared by single modes in non- symmetric Gaus- 
sian states with a well defined symplectic spectrum. We then 
perform an extensive numerical investigation which corrob- 
orates the measure of genuine multipartite entanglement in- 
troduced in Sec. |Vl and provides strong evidence that multi- 
partite entanglement is indeed strongly monogamous even in 
Gaussian states that are not symmetric under permutations of 
the modes. Conclusions and further outlook are discussed in 
Sec.lVnl 



II. PRELIMINARY DEFINITIONS 

A. Gaussian states and Gaussian operations 

Detailed accounts on the mathematical description of Gaus- 
sian states in the phase space formalism can be found in many 
references, e.g. |l9t l2l|]. In view of the subsequent analy- 
sis on quantum correlations, here it is sufficient to recall that 
Gaussian states of a bosonic A/'-mode CV system are com- 
pletely described in phase space (once the irrelevant first mo- 
ments are set to zero via local displacements) by the real, 
symmetric co variance matrix (CM) a of the second moments, 
whose entries are aij = l/2{{Xi^ Xj}) — {Xi){Xj). Here 
X = {xi,j9i, . . . , xat, Pat} is the vector of the field quadra- 
ture operators. The canonical commutation relations can be 
expressed in matrix form: [A^, Aj] = 2iftij, with the sym- 
plectic form Q = 0^=iCJ and 00 = — i^j = 1, 2. 
The CM a must fulfill the Robertson- Schrodinger uncertainty 
relation HI 

a^in>0 (1) 

in order to describe a physical state. Throughout the paper, 
a will be used indifferently to indicate the CM of a Gaussian 
state or the state itself. 

Unitary Gaussian operations U amount, in phase space, to 
symplectic transformations S. The latter are the transforma- 
tions that preserve the symplectic form, Q = S^ftS) acting 
"by congruence" on the CM (i.e. so that a ^ SaS^). Ac- 
cording to Williamson theorem f^T], the CM of a A^-mode 
Gaussian state can be always brought to diagonal form by 
means of a global symplectic transformation (this corresponds 
to the normal mode decomposition): 

W^aW^ = iy, (2) 

where Wo- G Sp{2N^ R) and v = diagji^fe, u^} is the 

CM of a fully uncorrected state corresponding to the tensor 
product of single-mode thermal states. The quantities {uk} 
are the so-called symplectic eigenvalues of the CM a. A pure 
Gaussian state is characterized by Vk = 1, V /c = 1 ... A", 
which implies det a = 1 and corresponds to the saturation of 
Ineq. ©• For mixed Gaussian states g with CM tr, the purity 
/i = Tr is computable as /i = 1 / Vdeta. 

For what concerns the characterization of bipartite entan- 
glement, the criterion of positivity under partial transposition 
(PPT) states that a Gaussian state a is separable (with respect 
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to a 1 1 bipartition) if and only if the partially transposed CM 
satisfies the uncertainty principle given by Eq. ([T]) [ 35, 36]. 
The suffix 'T" denotes the partial transposition, implemented 
by reversing the direction of time in the subspace of either one 
of the two subsystems in a composite, bipartite CV system 
1I35I1 . An ensuing computable measure of C V entanglement is 
then the logarithmic negativity [yf] Ej\f = In ||^^||i, where 
II • 111 denotes the trace norm. This measure is an upper bound 
to the distillable entanglement 1 38l ^|] of the state g. 

We employ a related but different measure of bipartite en- 
tanglement: the contangle 1 14], an entanglement monotone 
under Gaussian local operations and classical communication 
(Gaussian LOCC), that belongs to the family of 'Gaussian 
entanglement measures' HO, The principal motivation 
for this choice is that the contangle, unlike the logarithmic 
negativity, naturally enables a mathematical treatment of dis- 
tributed CV entanglement as emerging from the fundamental 
monogamy constraints (lEHsl]- The contangle Gr is defined 
for pure states as the square of the logarithmic negativity and 
it is extended to mixed states via the Gaussian convex roof 
construction |40l,l4p], that is as the infimum of the average 
pure-state entanglement over all decompositions of the mixed 
state in ensembles of pure Gaussian states. If ai\j is the CM 
of a (generally mixed) bipartite Gaussian state where subsys- 
tem i comprises one mode only, then the contangle Gr can be 
computed as IM liTIl 

^r(o-i|j) = Grio-'^^j) = g[d^j], g[x] = arcsinh^[Vx - 1], 

(3) 

where t7^|^.^ corresponds to a pure Gaussian state, and di\j = 

d{a^^j) = ydeto^ = y^deto^^, with tr^^^^^ the reduced 
CM of subsystem obtained tracing over the degrees of 
freedom of subsystem j (i). The CM a^^^^ denotes the pure 
bipartite Gaussian state which minimizes d{a^^.) among all 
pure-state CMs a^^. such that a^^. < a^j. If a^j is a pure 
state, then t7^|^.^ = ai\j, while for a mixed Gaussian state 
Eq. © is mathematically equivalent to constructing the Gaus- 
sian convex roof. For a separable state, d{a^^^^) = 1 which 
trivially means a vanishing contangle. The contangle Gr is 
completely equivalent to the Gaussian entanglement of forma- 
tion (4^, which quantifies the cost of creating a given mixed, 
entangled Gaussian state out of an ensemble of pure, entan- 
gled Gaussian states. Notice also that in general Gaussian 
entanglement measures and negativities are inequivalent, in 
that they may induce opposite orderings on the set of entan- 
gled, mixed non-symmetric two-mode Gaussian states HlW . 
This result reflects the fact that different measures of mixed- 
state entanglement are built on generally inequivalent opera- 
tional definitions and thus quantify different 'types' of entan- 
glement. 



B. Distributed entanglement and monogamy constraints 

Quantum entanglement, unlilce classical correlations, is 
said to be monogamous [12], meaning that it cannot be freely 



shared among multiple subsystems of a composite quantum 
system (see Ref. L42] for a brief topical review). In its orig- 
inal (we may call it "weak") form, the monogamy constraint 
imposes the following trade-off on bipartite entanglement dis- 
tributed among TV parties pi ... Pat, 

^pi\{p2...pn) y , (4) 

where £^ is a proper measure of bipartite entanglement (which 
may vary depending on the dimensionality of the system 
under consideration), and © is known as the generalized 
Coffman-Kundu-Wootters (CKW) inequality [13]. The left- 
hand side of inequality Q is the bipartite entanglement be- 
tween a probe subsystem pi and the remaining subsystems 
taken as a whole. The right-hand side is the total bipartite 
block entanglement between pi and each of the other subsys- 
tems pj^i in the respective reduced states. The non-negative 
difference between these two quantities, minimized over all 
choices of the probe subsystem, is referred to as the resid- 
ual multipartite entanglement: it quantifies the purely quan- 
tum correlations that are not encoded in pairwise form. In 
the simplest nontrivial instance of TV = 3, the residual en- 
tanglement has the meaning of the genuine tripartite entangle- 
ment shared by the three subsystems [13]. Such a quantity 
has been proven to be a tripartite entanglement monotone for 
pure three-qubit systems [^] , when bipartite entanglement is 
quantified by the tangle lfl3l] . and for pure three-mode Gaus- 
sian states, when bipartite entanglement is quantified by the 
contangle [14]. We recall that the general Ineq. (fflhas been 
established for spin chains (A/^-qubit svstems ) |[o[l5[l27l] and 
harmonic lattices (A^-mode Gaussian states) [14,76], with im- 
portant consequences for the structure of correlations of those 
many-body systems 0, |43l 113] . A proof of Ineq. for ar- 
bitrary quantum systems under arbitrary partitions has been 
also obtained llisi] when bipartite entanglement is quantified 
by the squashed entanglement [46], a measure which unfortu- 
nately lacks computability for almost all quantum states. 

III. STRONG MONOGAMY OF BIPARTITE AND 
GENUINE MULTIPARTITE ENTANGLEMENT 

In this Section we review the construction introduced in 
1I28I1 to refine the right-hand side of Ineq. © and to recur- 
sively extract a quantification of genuine multipartite entan- 
glement via a strong monogamy constraint on bipartite and 
multipartite quantum correlations (see Fig. \T\ for a pictorial 
representation). 

We investigate whether there exists a suitable generaliza- 
tion of the tripartite CKW-like analysis to an arbitrary num- 
ber N of modes, such that a genuine quantifier of TV-partite 
entanglement is naturally derived from a stronger monogamy 
inequality. We are motivated by the fact that the residual mul- 
tipartite entanglement emerging from the "weak" inequality 
^ includes all manifestations of i^T-partite entanglement, in- 
volving K subsystems at a time, with 2 < K < N. Hence, 
it severely overestimates the genuine A^-partite entanglement 
for N > 3. It is then compelling to introduce proper further 
decompositions of the residual entanglement, by performing 
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FIG. 1 : (Color online) Pictorial representation of a quantum system 
multi-partitioned into K subsystems {K — Q 'm the figure), each con- 
taining in general one or more elementary units (the total number of 
elementary units, depicted as spheres, is = 16 in the figure). The 
"weak" monogamy of bipartite quantum correlations imposes that 
the entanglement between a reference (yellow) subsystem and the 
rest of the system, exceeds the sum of the bipartite entanglements 
between the reference subsystem and each other (blue) individual 
subsystem. In Ref. |28] we have further decomposed this excess 
between the global block entanglement and the total bipartite entan- 
glement, into a sum of genuine multipartite entanglements involving 
3, 4, K parties at a time. The last contribution quantifies the truly 
K-partite quantum correlations shared among all parties. It is natu- 
ral to expect that this term is non-negative: such an argument yields 
a stronger monogamy constraint on the simultaneous distribution of 
both bipartite and multipartite entanglement. In Ref. |28] we proved 
the validity of such a decomposition and hence the existence of such 
a strong monogamy, for a system of N harmonic oscillators, one 
per party, in a globally permutation-invariant Gaussian state. Here 
we extend the strong monogamy framework to multi-partitions com- 
prising an arbitrary number of modes per party, retrieving a measure 
of genuine multipartite entanglement for larger classes of Gaussian 
states beyond symmetry constraints. 

finer investigations on the different ways in which a subsys- 
tem can be entangled with the entire group of the remaining 
N — 1 subsystems. 

Quite naturally, one can argue that a single subsystem 
can share individual pairwise entanglement with each of the 
other subsystems; and/or genuine three-partite entanglement 
involving any two of the remaining subsystems, and so on, 
until it can arrive to be genuinely A/^-party entangled with all 
the remaining subsystems. We thus introduce the hypothesis 
that these contributions are well defined and mutually inde- 
pendent, and verify a posteriori that this is indeed the case. 
Namely, we wish to verify whether in multipartite states, en- 
tanglement is strongly monogamous in the sense that the fol- 
lowing equality holds |28]: 

+ ... + £;pib2|...b7v ^ (^5) 

where E^^^p^ is the bipartite entanglement between parties 
1 and j, while all the other terms are multipartite entangle- 
ments involving three or more parties. The last contribu- 
tion in Eq. ^ is defined implicitly by difference and repre- 



sents the residual A^-partite entanglement. For generally non- 
symmetric states, one can define a set of residual entangle- 
ments £;^l^*2 1 • • • AT which may possibly differ quantitatively 
depending on the probe system with respect to which the 
total one-to-A^ bipartite entanglement is decomposed. In or- 
der not to overestimate the entanglement shared by all the in- 
dividual parties, it is thus fairly natural to define the genuine 
residual A/^-partite entanglement as the minimum over all the 
permutations of the subsystem indexes, 

^Pl\P2\---\PN = j^^j^ ^Pil\Pi2\---\PiN ^ (^5) 

{ii,...,iAr} 

We anticipate that such a definition is also necessary as, for in- 
stance, the average or the maximum of the residual entangle- 
ments fail to satisfy the strong monogamy condition. For fully 
symmetric states, there is no need to perform a minimization 
as the residual entanglements are all equal and independent on 
the probe subsystems. 

All the contributions to multipartite entanglement appear- 
ing in Eq. ([5]) (except the last one) involve K parties, with 
K < N, and are of the form ^^^i b2 1 • • • bi^ . Each of these terms 
is defined by Eq. ([5]) when the left-hand- side is the 1 1 (iC — 1) 
bipartite entanglement £;^^iI(^^2 . pk) A^-partite entangle- 
ment is thus, at least in principle, computable in terms of the 
known i^T-partite contributions, once Eq. ([5]) is applied recur- 
sively for all = 2, . . . , A/" - 1. To assess ^^^ib2|...bAr a 
proper quantifier of A/^-partite entanglement, one needs first to 
show its non-negativity on all quantum states. This property in 
turn implies that Eq. ([5]) can be recast as a sharper monogamy 
inequality 

^P1\P2\...\PN > 0^ (7) 

constraining both bipartite and genuine K-partite (K < N) 
entanglements in A/'-partite systems. Such a constraint on dis- 
tributed entanglement thus represents a strong generalization 
of the original CKW inequality Q [13], implying it, and re- 
ducing to it in the special case A^ = 3. 

In Ref. [28] we have further simplified the construction of 
Eq. dl]) under the assumption ("full symmetry") that the states 
under consideration are invariant under any permutation of the 
subsystems. Symmetric states play in fact a prominent role in 
multiparty quantum information applications ||47|] and theo- 
retical entanglement investigations |48i|^. For permutation- 
invariant Gaussian states we obtained a closed expression for 
the genuine multipartite entanglement distributed among sin- 
gle modes in terms of an alternating sum of bipartite block 
entanglements only [28]. We were moreover able to prove 
analytically that multipartite entanglement distributed among 
single-mode parties is strongly monogamous in arbitrary A^- 
mode Gaussian states that are fully symmetric under arbitrary 
permutations of the modes ll28ll . We eventually extracted a 
proper, computable measure of the genuine multipartite en- 
tanglement shared by the N single modes in fully symmetric 
Gaussian states. We proved that the A^-partite residual en- 
tanglement is monotone in the optimal nonclassical fidelity 
of A/^-party teleportation networks with symmetric Gaussian 
resources 1101123], and thus acquires an operational interpre- 
tation and a direct experimental accessibility ||28|] . 




FIG. 2: (Color online) Pictorial representation of the process of unitary localization and delocalization of multipartite entanglement in multi- 
symmetric Gaussian states. The left picture represents a (pure or mixed) Gaussian state of N modes (N = 16 in the figure), all of them in 
general correlated, described by a covariance matrix a. The A^-mode system is multi-partitioned into K < N "molecules" Mi, ... Mr 
(i^ = 6 in the figure). Each molecule Aij consists in a block of mj modes, each of them having the same reduced covariance matrix (the 
same color in the figure). By definition, multi- symmetric states are invariant under mode permutations within a single molecules. There exists 
a symplectic transformation Sioc = Wi W2 ... Wk, acting locally on each molecule, which transforms the original state into a new 
one, described by a covariance matrix a' = SiocCrSfoc^ such that all correlations are concentrated among K modes only ("nuclei"), one 
per molecule, as illustrated in the right picture. Precisely, cr' corresponds to a generally entangled state of the K nuclei, tensor {N — K) 
single-mode thermal states which do not bear any correlation with each other and with the nuclei. Any bipartite and multipartite form of block 
entanglement shared by the molecules in the state cr, is then converted {localized) into an effective entanglement among single modes (realizing 
the transition from left to right in the figure). The localization is reversible, as Sioc amounts to a local unitary operation at the level of Hilbert 
space, hence the entanglement can be redistributed (delocalized) among all modes again by applying the inverse symplectic transformation 
'^lol ^ (realizing the transition from right to left in the figure). 



Here we discuss a more general case in which fully sym- 
metric A^-mode Gaussian states are partitioned into K blocks 
of modes — "molecules" — of arbitrary size. We then ana- 
lyze the validity of the strong monogamy and the quantifica- 
tion of the genuine multipartite entanglement among the K 
molecules. A key result on which our analysis is built con- 
cerns the possibility of transforming by local (i.e. acting on 
each molecule) unitary operations, such a i^T-partite block en- 
tanglement into the i^T-partite entanglement of an equivalent, 
generally non- symmetric, Gaussian state of K modes only. 
We will now show that this property holds indeed and is a 
consequence of a general result in symplectic analysis. 



sentation). With this notation we refer to a Gaussian state of 
a A/^-mode CV system, multi-partitioned into K subsystems 
{K < N). Each subsystem Mj (j = 1, . . . , K) is termed a 
"molecule" and consists in a block of rrij modes. The multi- 
symmetry (specifically, i^^-symmetry) requirement means that 
the state is invariant under the permutation of any two modes 
within each molecule. This implies, in particular, that the 
rrij modes comprising the j-th molecule all have the same 
marginal covariance matrix. Explicitly, 



IV. UNITARY LOCALIZATION OF MULTIPARTITE 
GAUSSIAN ENTANGLEMENT 

Let ^''^j^^ be the CM of a multi-symmetric, gen- 

erally mixed Gaussian state (see Fig. [2] for a pictorial repre- 
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Here 5]j is the reduced (fully symmetric) CM of the block 
of modes forming the molecule Mj, while the j matrices 
encode the correlations between two molecules Mi and Mj. 
Each lower-case Greek letter corresponds instead to a two by 
two matrix. 

In Ref. [32] we have investigated the particular instance of 
bi-symmetric, bipartite Gaussian states (i.e. the case K = 2). 
There, we found that nonlocal correlations in such states are 
unitarily localizahle: There exist local (i.e. acting indepen- 
dently on each party) and unitary (i.e. symplectic at the CM 
level) operations that transform the state into a correlated state 
of two modes only (one per party), tensored with — 2 uncor- 
rected single-mode states. This result stems from the degen- 
eracy of the symplectic spectrum of the reduced (fully sym- 
metric) state of each party |l32t |48D. Indeed, the proof im- 
mediately extends to the more general, multipartite case that 
we investigate in the present work. The crucial observation 
is that the local unitary operations enabling the localization 
are exactly the Williamson transformations Wj which realize 
the normal mode decomposition, Eq. Q, in each molecule 
M.j. Therefore, they only depend on the molecule on which 



they act. This entails that localizing operations Wj applied 
to different molecules commute, and all molecules can be si- 
multaneously diagonalized by a "local" (with respect to the 
given K-partition) symplectic operation 5'/oc ^ 5'p(2A/', R), 
obtained as the direct sum of the individual Williamson oper- 
ations Wj\ Sloe H^i © V^2 © . . . © Wk- Trivially, each 
two-molecule reduced state of a multi- symmetric state of the 
form Eq. ([5]) is a bi-symmetric state that can be immediately 
and unitarily localized by the action of Sioc- Specifically, one 
then has that the original CM cr^J^]^ _ where in prin- 
ciple all modes could have been correlated with each other, is 
transformed into a new CM 



K:{mi,...,mK} 



The new CM cr'^) . corresponds to a Gaussian state 

in which all correlations are concentrated (localized) among 
K modes only — that we term "nuclei" — one per molecule, 
and the additional (N—K) modes are tensored out. Explicitly, 
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where we have denoted by a the effective K-modt CM of the 
nuclei, while the CMs Uj represent the uncorrelated single- 
mode thermal states. The elements of the CM a' can be de- 
termined in general from the local (molecular) symplectic in- 
variants of a [|32l1 . 

^From the previous analysis we see that the property of 
permutation invariance within each block of modes implies a 
powerful transformation of multipartite Gaussian states, such 
that the correlations (including entanglement) shared by the 
multimode parties can be computed equivalently as correla- 
tions among single modes in an effective two-mode Gaussian 
state. This localization is perfectly reversible, being enabled 
by the unitary (symplectic) local operation Sioc, such that the 
action of 6'^"^ on a' just delocalizes correlations and redis- 
tributes them among the various modes of each molecule, re- 
constructing the original multipartite Gaussian state a. 

An important consequence of the property of unitary local- 
ization is that for the relevant subset of Gaussian states that 
obey the constraint of multi- symmetry, the computation of the 
multipartite molecular entanglement is greatly simplified by 
reducing it to the entanglement shared by single modes. For 
general, non-symmetric multipartite Gaussian states, the min- 
imal form of the CM that can be obtained by local unitary 



operations involves instead correlations among more than one 
mode per molecule. Specifically, by virtue of the general de- 
composition proven in [50], there remain irreducible mode- 
wise and pairwise correlations among the elements of each 
molecule^. It is known that the pairwise correlations vanish 
for pure, generally non-symmetric Gaussian states jsi] ]. Here, 
we have further proved that they vanish as well for generally 
mixed, multi- symmetric Gaussian states, and that the mode- 
wise correlations only involve one nuclear mode per party. 

We remark that thanks to the unitary localization of Gaus- 
sian entanglement, already restricted to the case of bipartite 
bi-symmetric states fs^], one can immediately conclude that 



By mode-wise correlations we mean those contained in a state in which 
single modes each pertaining to a different molecule are correlated. There 
may be in general a collection of such mode- wise correlated states, in- 
volving independently the modes belonging to each molecule. By pair- 
wise correlations one means instead those correlations contained in a state 
in which pairs of modes belonging to the same molecule (being them- 
selves correlated) are inter-correlated with pairs of modes belonging to 
other molecules. More than one pair of modes per molecule can exhibit 
such correlations in general states. 
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entanglement distillation with correlated copies (S^ is pos- 
sible with Gaussian states and Gaussian operations. This is 
in sharp contrast to what happens with the standard paradigm 
of entanglement distillation, where it is known that from sev- 
eral independent and identically distributed copies of mixed 
bipartite Gaussian states one cannot obtain less copies of more 
entangled two-mode squeezed states by performing Gaussian 
LOCCs alone [53]. We have just shown instead that the pres- 
ence of correlations between the copies (testified by a global 
fully symmetric pure state of the form Eq. ^ with nonzero 
Tij's) enables Gaussian entanglement distillation by simple 
unitary Gaussian local operations. Finally, it is worth men- 
tioning that results similar to the ones presented in this Section 
have been demonstrated concerning the possibility of extract- 
ing GHZ-type (GHZ stands for Greenberger-Horne-Zeilinger 
[54]) entanglement by local unitary operations performed on 
multipartite A^-qubit stabilizer states ll55h . which are some- 
how regarded as the discrete variable analogues of Gaussian 
states 15 6] . 



V. GENUINE MULTIPARTITE BLOCK ENTANGLEMENT 
OF FULLY SYMMETRIC GAUSSIAN STATES: 
MOLECULAR HIERARCHY 

Fully symmetric Gaussian states are multi- symmetric with 
respect to any global multi-partition of the modes. The CM of 
such states can be written in standard form as lISj, Hi 



rW = 

fs 



/I3 C 
C /3 

'■■ c 
VC ••• 



(10) 



with (3 = diag(6, 6) (b > 1) and = diag (^i, 22)- 
Pure fully symmetric Gaussian states of N modes 
are characterized by zi = [(6^ — 1)(A'" — 2) + 
V^(62 - l)(62jV2 - (iv - 2y)]/[2b{N - 1)], Z2 = [{P - 
l){b^N^ - (AT - 2)2)]/[26(iV - 1)]. 



1){N - 2) + 

A Mathematica code instruction to construct their CM for 
any number N of modes is provided by Eq. (A. 2) in the 
Appendix. 

Fully symmetric Gaussian states, Eq. ([TOb , are key re- 
sources for essentially all the multiparty CV quantum infor- 
mation protocols implemented so far [20], including the tele- 
portation network llslll^], and the Byzantine agreement |[3QIi . 
They can be experimentally prepared by sending a single- 
mode squeezed state with squeezing in momentum and 
(A^ — 1) single-mode squeezed states with squeezing Tp in 
position, through a network of A" — 1 beam-splitters with 
tuned transmittivities, as shown in detail in Refs. 
Up to local unitaries, the CM of Eq. ([TOb can be parameter- 
ized equivalently only by considering the average squeezing 
r = (r^n -\- rp)/2 needed in the preparation of the state. In 
fact [code: Eq. (A.3)], 

2 _ (A- - 2)A^ + 2{N - 1) cosh(4r) + 2 



Recall that for 6 > 1, i.e. \r\ > 0, the associated fully sym- 
metric Gaussian state is fully inseparable. In general, the de- 
terminant of the reduced L-mode CM cr^^-' of a fully symmet- 
ric A'-mode pure Gaussian state is given by [code: Eq. (A.4)] 



dettr^^^ = [2L^ -2NL^2{N-L)cosh{4r)L^N^]/N'^. 

(11) 

Here, the CM tr^^^ denotes a fully symmetric Gaussian state, 
which is mixed as a consequence of the operation of trac- 
ing over N — L modes. Recalling that the determinant of 
the CM a is related to the purity /i of a Gaussian state by 
/i = (det tr)"^/^, we see that the purity of a^^-' strictly de- 
creases with increasing A", and is concave as a function of L 
with a maximum at L = A'/2 (A^ even) or L = (A' ± l)/2 
(N odd). 

/,From now on, by ^''^j^^ we will denote a pure, 

fully symmetric Gaussian state of N modes, multi-partitioned 
into K molecules, respectively made of mi, ... , rrik modes. 
Here we allow rrij = M to be in principle smaller than 
A", thus considering in general a mixed, symmetric multipar- 
tite entangled state of the K molecules. The genuine multi- 
partite entanglement, quantified by the residual contangle as 
defined via the strong monogamy constraint Eq. Q, has been 
computed in the case of rrij = 1 Vj | 28], i.e. of mono-modal 
molecules ("atoms"). We will now show how to evaluate the 
multipartite entanglement shared by molecules of arbitrary 
size (multi-modal molecules). 

Recursive applications of Eq. ^ eventually reduce the task 
to the computation of all the bipartite mi\mj entanglements. 
As we have shown in Sec. [iVl the bipartite entanglement be- 
tween rrii modes and rnj modes in a fully symmetric Gaussian 
state can be concentrated by local unitary operations onto two 
modes only (unitary localization) [32]. This effective two- 
mode Gaussian state turns out to be a mixed state of partial 
minimum uncertainty, the so-called GLEMS (Gaussian Least 
Entangled Mixed State) | 57]. A two-mode GLEMS is defined 
by having the smallest symplectic eigenvalue equal to 1. This 
means that one of its normal modes is in the vacuum state 
(which is trivially a Heisenberg minimum-uncertainty state) 
and all the mixedness is concentrated in the remaining normal 
mode. Moreover, a two-mode GLEMS is completely speci- 
fied by its global (two-mode) purity and its two local (single- 
mode) purities, respectively given by 



l/c 
1/a 
1/6 



(12) 



where the determinants are defined by Eq. ( fTTI) . The en- 
tire family of Gaussian measures of entanglement, including 
the contangle Gr L14.] . can be computed in closed form for 



10 



GLEMS m, yielding [with the notation of Eq. ©] 



c < ci; 



(C2-1)^ ' 

.1, c > C2; 

(13) 

H ere ci = [(a' - + 2(a2 + b^) + 

v/(a2 -62)2 + 8(a^ + 62)|a' - 62|]V2/[2(a2 + b% 
C2 = + 6^ — 1, and S = {a — b — c — l){a — b — c -\- 
l)(a + 6- c-l)(a + 6- c+l)(a-6 + c-l)(a-6 + c + 
l)(a + 6 + c - l)(a + 6 + c + 1). Eq. ^ is coded in the 
Appendix as Eq. (A.5). 

Therefore, the bipartite contangle between any two 
molecules Mi and Mj which are, respectively, -modes 
and mj -modes sub-blocks of a fully symmetric pure A^-mode 
Gaussian state, is analytically determined by the expression: 



rrii \ mj 



(14) 



arcsinh^ 




det a^^} , 



det a^^} , 



det a 



(AT) 



together with Eqs. ([TTHT31) . 

We are now finally in the position to evaluate Eq. Q 
and determine the genuine iiT-partite entanglement 
Qres{^W \ (iistributed among K molecules in 

T \ mi |m2 1 ..-I'mx ^ ^ 

a globally pure or mixed fully symmetric Gaussian state of 



K 



rrin 



M < N modes. Identifying each subsystem 



with a molecule, and adopting the contangle as a measure of 
bipartite entanglement, Eq. Q can be rewritten as 



r V mi|m2|...|mx^ 

= G.(a(^) ) 

^ V mi > . . , m.,. / 



(15) 



J2>1 

EQresf (N) 
r \ mi\mj^\mj. 

j3>j2>l 



E 



jK-l>jK-2>--->j2>l 



(AT) 

mi|mj2 I- 



.\mj 



Here we have relabeled the molecules so that their sizes are 
sorted in ascending order, mi < 1712 < ... < tuk, and we 
have used the fact that the multipartite entanglement is mini- 
mized when the probe subsystem has the smallest determinant 
among the involved parties 1 14]. Eq. ilSi can be easily solved 
by using a computer language that supports recursion. For 
this purpose, as well as for other computational tasks of this 
Section, we have written a Mathematica instruction program 
[Eq. (A. 9) in the Appendix] that evaluates exactly and analyt- 
ically Eq. ([T5]) for any fixed value of the squeezing r, of the 
number of modes N, of the number of partitions K, and for 
each given partition {mi , . . . , tuk}. 



The resulting formalism is quite compact, notwithstanding 
that one is dealing with a complicated combination of nested 
contributions to multipartite entanglement with alternating 
signs, and it is by no means apparent a priori, even though 
intuitively expected, that the quantity G^^*(o-^j^^| |^^) 
is non-negative and thus verifies the postulated property of 
strong monogamy of multipartite entanglement, Ineq. (|7]). Al- 
though achieving a full analytical proof in the most general 
situation appears to be a formidable task, we can test the va- 
lidity of our construction in a wide variety of different situa- 
tions. In this way, we will single out the peculiar features of 
multipartite Gaussian entanglement, and its scaling with the 
number and/or the size of the molecules as well as with the 
global size of the system and the number of modes that are 
traced out. As we will show, all investigations confirm the 
non-negativity of the proposed measure of genuine multipar- 
tite entanglement and hence support the validity of the strong 
monogamy decomposition. In the following, we present a sys- 
tematic discussion of the various cases and situations consid- 
ered. 



A. Atomic multipartite entanglement 

By the term "atomic entanglement" we denote the one 
shared by individual modes, i.e. by molecules with rrij = 1 
Vj = 1, . . . , iC. In this case, one trivially has rrij = 

M = K, and the associated K-mode state tr^j-^ -^^ is 
obtained from a pure, fully symmetric Gaussian state of N 
modes, Eq. (fTOb , by tracing over N — K modes. In this spe- 
cific case, one can derive an explicit formula for G^^^ and its 
non-negativity can be proven analytically for any N, K, and 
r (the single-mode squeezing) [ 28]. We review this important 
result first, skipping the technical details, for the sake of com- 
pleteness. The genuine iiT-partite entanglement reads [code: 
Eq. (A.7)] 




(16) 



X arcsinh 



2Vi^ - 1 - jsinh(2r) 



and one can show that G^^^(tT^j-^ ^y) > (see Fig. O, 
proving the strong monogamy Eq. dJj) of bipartite and multi- 
partite entanglement distributed among_single modes in a fully 
symmetric Gaussian state ll28h . Ref. 11281] provides as well an 
operational interpretation for G^^^ in the case of pure states 
(K = N) by showing its monotonic dependence on the op- 
timal fidelity of CV teleportation networks ||T3, [H] for any 
N. Within the scope of this paper, we are interested in study- 
ing the scaling of multipartite entanglement with the relevant 
physical parameters characterizing the multimode Gaussian 
states. 
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FIG. 3: (Color online) Genuine multipartite entanglement G^^-^^ of 
pure, fully symmetric A^-mode Gaussian states cr^j^ plotted 

as a function of N (up to 10^) and of the average squeezing r needed 
in the preparation of the state. All the quantities plotted are dimen- 
sionless. 



^From Fig. [3] we observe that G^^* increases monotoni- 
cally with the average squeezing r, while it decreases with the 
number of parties K, eventually vanishing in the field limit 
K ^ oo. For mixed states {K strictly smaller than N), G^^* 
decreases with the number N — K of the modes that have 
been traced out, i.e. with the mixedness. The monotonically 
decreasing dependence on K can be understood as follows. 
According to the strong monogamy decomposition Eq. (|5]), as 
the number of modes increases, the residual non-pairwise en- 
tanglement can be encoded in so many different multipartite 
forms, that the genuine K-partite contribution is strongly sup- 
pressed (frustrated). 



B. Molecular multipartite entanglement versus squeezing: 
promiscuous entanglement distribution 

^From now on we consider general Gaussian states 
(T^J^]_ _ Recall that N is the total number of modes 

K:|mi,.. .,771x1 

for a fully symmetric pure Gaussian state cr^^J , Eq. ([TOb , from 



which a 



(AT) 



M 



F r .is obtained by tracing over N 

modes (here M = rrij). The K molecules, i.e. the parties 
whose shared genuine multipartite entanglement G^^^ is being 
computed via Eq. ([15]), are formed by mi , ... , ttik modes, re- 
spectively. 

A unifying trait in our investigations is that G^^^ is found 
to be non-negative on any partition (for a comparison with ex- 
tended numerical tests, see Sec.lVll). Moreover, it is a mono- 
tonically increasing function of the squeezing r that enters in 

the engineering of the state a^j^-' . However, there is a signif- 
icant difference between the case in which the state of the K 
molecules is globally pure (M = N) and the one in which it 
is mixed (M < N). In the first instance, for any N, K, and 
individual molecule sizes m^'s, G^^^ increases unboundedly 
with r and diverges as r ^ oo (which is accompanied by a 
divergent mean energy). For a mixed molecular state, instead. 



for any number and size of the molecules the genuine multi- 
partite entanglement remains always finite and saturates at a 
constant asymptotic value in the limit of infinite squeezing. 
Examples of these two different behaviors will be illustrated 
and discussed in the subsequent graphics. 

Extensive numerical analysis shows that, as soon as r > 0, 
all possible multi-partitions ^''^j^^ exhibit genuine 

multipartite correlations, even though the presence of mixed- 
ness limits the maximal amount of achievable entanglement. 
This coexistence of bipartite, tripartite, . . . , genuine i^-partite 
(etc.) entanglement, together with the fact that all the var- 
ious contributions are monotonically increasing functions of 
each other (being all increasing as a function of r), defines 
the so-called "promiscuous'' structure of entanglement shar- 
ing in Gaussian states Promiscuity is a general feature 
of distributed entanglement that arises with increasing dimen- 
sion of the Hilbert space llssll . It is fully displayed in Gaussian 
states of CV systems where despite the firm holding of strong 
monogamy constraints, there remains a certain freedom in the 
way quantum information can be shared. 

^From an operational point of view, we can conclude this 
subsection by stating that multi-partitioning the global N- 
mode Gaussian system without discarding any mode is the 
best choice in order to take full advantage of the squeezing re- 
source and to exploit the infinite entanglement capacity which 
characterizes CV systems. We will now investigate the de- 
pendence of multipartite entanglement on the specific way in 
which a CV system is multi-partitioned. 



C. Molecular multipartite entanglement vs. number of modes 



Here we analyze how the iiT-partite molecular entangle- 



ment in the state o-^^) . is affected by a variation of 

N. We have to distinguish between two possibilities: (i) a 
variation of N with an accompanying variation of the m^'s, 
i.e. a global re-scaling of the system; (ii) a variation of N 
keeping fixed molecular sizes, i.e., a variation of the mixed- 
ness. 

(i) In this case we observe the remarkable phenomenon that 
the genuine multipartite molecular entanglement of symmet- 
ric Gaussian states is scale invariant. Let us consider the state 
o- . , and construct another state cr 1^^^ . 

K:{mi,...,mK } K:{s mi . .,s ttlk } 

where both the total number of modes N and the sizes of the 
individual molecules mi , . . . , tuk (and hence their sum M) 
are re-scaled by an arbitrary integer factor s. It is immediate 
to see that the genuine multipartite entanglement G^^^ among 
the K molecules is equal in both states. This follows by recall- 
ing from Eqs. ([T2l - [T4b that G^^^ depends ultimately only on 

the reduced determinants det tr^^^ of L-mode sub-blocks of 

the A/'-mode state, and by observing that det cr[^^-' in Eq. ([TT]) 
is independent of the scale factor s. Therefore we see that size 
does not matter for what concerns multipartite Gaussian en- 
tanglement, if the partition structure is preserved. This is rel- 
evant in view of the practical exploitation of Gaussian states 
for communication tasks |20]: Adding redundance, e.g. by 
doubling the size of the individual subsystems, yields no ad- 
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FIG. 4: (Color online) Genuine multipartite entanglement G^^^^ in 
the state 2 3 4} ^ function of the squeezing r. The plot re- 
ports the entanglement shared by four molecules made respectively 
of 1, 2, 3, and 4 modes in 10-mode fully symmetric Gaussian states, 
obtained, respectively, as reductions of a pure, fully symmetric A^- 
mode Gaussian state with = 10 (solid red line, no modes are 
traced out), = 11 (dashed magenta line), = 12 (dot-dashed 
green line), A^ = 13 (dot-dot-dashed cyan line), and A^ = 14 (dot- 
ted blue line). The genuine multipartite molecular entanglement de- 
creases with increasing N, i.e. with the number of discarded modes. 
In the pure- state instance (no modes being discarded) it diverges as 
r ^ 00. All the quantities plotted are dimensionless. 



vantage for the multiparty-entangled resource. 

It is interesting to observe that an analogous property of 
scale invariance holds as well for the set of hierarchical, multi- 
component measures of geometric entanglement in arbitrary, 
fully symmetric pure states of A^-qubit systems JH]. The fact 
that two different measures of genuine multipartite entangle- 
ment evaluated on two very different systems (A/^-qubit sys- 
tems and A/'-mode CV systems) exhibit the same structure of 
scale invariance hints at the intriguing possibility that multi- 
partite entanglement is indeed scale invariant on all fully sym- 
metric states of arbitrary multi-party quantum systems. 

(ii) On the other hand, if the number of modes N is 
increased to A^' > A^ without re- scaling the size of the 

molecules, this results in a state cr^!^ _ . which is al- 

' K :|mi ,...,m_fc| 

ways less pure than the state cr^j^^ mx}' ^^^^ because 
we are increasing the number of modes which are traced out. 
As it could be intuitively expected, in this case the genuine K- 
partite molecular entanglement decreases monotonically with 
N' . An example of this behavior when scaling the number of 
modes is illustrated in Fig.lH 



D. Molecular multipartite entanglement vs. number of 
molecules 




r 



0.2 0.4 0.6 0.8 1 

FIG. 5: (Color online) Molecular entanglement hierarchy in a pure 
fully symmetric M-mode Gaussian state with M = 12. From left 
to right: the bipartite entanglement between two 6-mode molecules 
(yellow line) is strictly stronger than the genuine tripartite entangle- 
ment between three 4-mode molecules (magenta line), and so on, 
while the atomic genuine 12-partite entanglement (red line) is the 
weakest one. The bipartite entanglement between any two modes in 
the reduced mixed state obtained by tracing over 10 modes is plotted 
as well as reference (dashed black line). All the quantities plotted are 
dimensionless. 



with M large enough to allow the system to be partitioned 
into Ki ox K2 • • • or Km equally sized molecules^ as well 
as into M atoms (each molecule formed by one mode). The 
integer M can be chosen as a common multiple of the set 
{i^Ti, . . . , Km}, in such a way that in each instance the sys- 
tem is partitioned into Ki molecules made of M/Ki modes. 
This is an interesting case to study as we may imagine the ex- 
perimental situation in which a given multimode symmetric 
Gaussian state has been prepared, and an option is available 
whether to split and distribute it to few or to many parties, in 
order for them to share as much entanglement as possible. 

By virtue of the scale invariance of G^^* (see Sec. IV Cb , 
the genuine multipartite entanglement shared by Ki even- 
sized molecules in the state o-^xJiM/Ki M/Ki} ^^^^ 
the genuine multipartite entanglement shared by Ki modes 
in the state ^^^^ji ^j, i.e. to the atomic entanglement of 
the latter state. We have shown in Sec. IV Al that the atomic 
entanglement, Eq. ([T6b , monotonically decreases with an in- 
creasing number of entangled modes, Ki in this case. These 
two observations allow to conclude that the genuine multipar- 
tite molecular entanglement in the state o-^x\m/k- m/k } 
is monotonically decreasing with an increasing number Ki 
of molecules. As a general rule, therefore, a smaller number 
of larger molecules shares strictly more entanglement than a 
larger number of smaller molecules: The maximum entangle- 
ment is always obtained by bisecting symmetrically the sys- 
tem and distributing it to two parties only. An example of this 
hierarchy of molecular entanglements is illustrated and dis- 



Next, we keep both N and M fixed, and study the depen- 
dence of G^^^ on the number K of molecules in which the 
M modes in the state ^''^j^^ are partitioned. To vi- 

sualize the problem, suppose to have pure states (M = A^) 



^ Although our findings do not rely on these assumptions, we consider pure 
states and molecules with the same size for simplicity in order to enable 
a more direct comparison, since the multipartite entanglement is further 
affected by the degree of mixedness (see previous subsections) and by the 
size of the individual molecules (see next subsection). 
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cussed in Fig.[5l 



E. Multipartite entanglement vs. size of the molecules 

Finally, we keep N, M, and K fixed. We are thus deal- 
ing with a well defined fully symmetric Gaussian state of M 
modes, with an assigned degree of mixedness determined by 
the number A^ — M of discarded modes, and, importantly, with 
a well defined number of molecules K in which the M modes 
are partitioned. This instance corresponds to the experimen- 
tal situation in which a protocol involving K parties has to 

be implemented, and cr^x-{mi ttik} resource state that 
needs to be employed. There is a freedom left in the choice 
of the individual sizes ruj of each molecule (constrained to 
rrij = M), which corresponds to the freedom in assign- 
ing a certain number of modes to each party. The distribution 
can easily break the symmetry of the original state whenever 
the different parties are assigned a different number of modes 
rrij . We thus wish to investigate how the genuine multipartite 
entanglement shared by the K molecules, Eq. ([15]), depends 
on the partition {mi, . . . , m^}, i.e. on the specific molecular 
size. 

Determining all the possible independent partitions of M 
modes into K subsets is a combinatorial problem. Here only 
the individual values of the integers rrij are relevant, not their 
order (e.g. {1,2,3} and {3, 2, 1} count as the same partition). 
It is thus convenient to introduce a ranking of the different 
partitions in order to be able to write a sorted list, without rep- 
etitions, of all the possible sets {mi, . . . , m^} given any inte- 
ger value of M and K. We denote by Vm,k = {Pm,k}^=i'' 
the list of partitions, whose i-th element is of the form 

P^M K ~ {^1^5 ' ' ' •> ^kI' cardinality of Vm,k is the 
integer Hm,k- The ranking we choose is very natural. Let 

Vm k the integer decimal number whose digits are given 

by the m^^^'s in base M, i.e. rj^M k ~ Ylf=i'^^j^ • 

Then the elements ^ of the list Vm,k are sorted in such 

a way that the corresponding ranking numbers r^y ^ are 
in ascending order. A simple way to construct the sorted 
list of partitions of M modes into K molecules is given 
by our Mathematica instruction KSortedPartitions[M , K], 
coded as Eq. (A. 8) in the Appendix. For instance, Vs 4 = 
{{1,1, 1,5}, {1,1, 2,4}, {1,1,3, 3}, {1,2, 2, 3}, {2, 2,2; 2}}. 

Having determined the combinatorial structure of the possi- 
ble independent partitions, we can investigate the dependence 
of the genuine multipartite entanglement on the molecular 
partitioning for a rich variety of pure and mixed fully sym- 
metric Gaussian states with different assigned values of TV, 
M, and K. Some representative instances are illustrated in 
Fig.[6l After an extensive numerical analysis, we have found, 
quite interestingly, that G^^^ appears to increase monotoni- 
cally with increasing ranking of the sorted list of all the pos- 
sible partitions. Specifically, the most inefficient way to par- 
tition the system is always the most unbalanced one that as- 
signs a single mode to all parties but one, and all the remaining 
modes to the last party. On the contrary, the best strategy to 



distribute the maximum if -partite entanglement corresponds 
to a balanced partitioning, in which all the parties have almost 
the same number of modes (i.e. all equal m^'s or the closest 
possible setting, depending on M being even or odd and on 
its ratio with K). Between these two extremal cases lie all 
the intermediate situations in which, gradually, the molecules 
tend to assume similar sizes (following the ranking of the list 
Vm,k), and the genuine iC-partite entanglement increases ac- 
cordingly. 

The overall picture shows that G^^^ exhibits a hierarchi- 
cal structure even with respect to the way the system is par- 
titioned. In fact, such a picture matches the intuition that 
molecules with similar sizes yield a stronger entanglement. 
It is quite natural indeed to expect that if the different par- 
ties share high genuine multipartite entanglement, they should 
carry rather similar amounts of quantum information (quan- 
tifiable by the local purity), which is precisely the case for 
molecular partitions with high ranking in our sorted list. 



VI. NUMERICAL TEST OF STRONG MONOGAMY. 

GENUINE MULTIPARTITE ENTANGLEMENT OF 
EFFECTIVE NON-SYMMETRIC GAUSSIAN STATES 

Throughout the previous Section we investigated the gen- 
uine multipartite entanglement G^^^ in states of the form 

^^K-{mi ttik}' entanglement shared by K molecules 

which are sub-blocks of a fully symmetric, pure or mixed 
Gaussian state. Although the sizes of the molecules were in 
general different, the severe symmetry constraint on the orig- 
inal state might induce to conclude that the strong monogamy 
of Gaussian entanglement is a feature limited to permutation- 
invariant states^. Actually, the unitary localizability of mul- 
tipartite entanglement (see Sec. HVb , is not only a key ingre- 
dient for the determination of G^^^ via Eq. ([15]), but also en- 
ables an immediate re-interpretation of the previous calcula- 
tions in terms of the genuine multipartite entanglement shared 
by K single modes (nuclei) in effective non- symmetric pure 
or mixed i^-mode Gaussian state <j x:{mi,...,mK}'' 

Eq. ©. 

In particular, in our investigations cr is the CM of fully 
symmetric Gaussian states, whose symplectic properties are 
known [ 32, 48]. As a consequence, the nuclear state a shares 
the distinctive property that all of its symplectic eigenvalues 
but one are equal to 1 and the non-unit symplectic eigenvalue 
is trivially equal to V det a, i.e. to the inverse of the purity 
of the state. States of this form can be regarded as partially 
saturating the Heisenberg uncertainty relations Ineq. ([B and 
are thus the direct multipartite generalization of the two-mode 
GLEMS [57] already introduced in Sec.|Vl We will therefore 
refer to state d-^J^\ . as a "multi-GLEMS" or K-mode 

K :|mi ,...,mx| 

GLEMS. 



^ Notice however that a family of non- symmetric four-mode Gaussian states, 
not related to the classes of states considered here and exhibiting an unlim- 
ited entanglement promiscuity, has also been proven to obey the strong 
monogamy constraint |58J. 



14 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.0 0.2 04 0.6 0.8 1.0 1.2 

r r 

(c) (d) 



FIG. 6: (Color online) Plot, as a function of the squeezing r, of the genuine iT-partite entanglement G^^^ in A^-mode pure fully symmetric 
Gaussian states, multipartitioned into K molecules of varying sizes {mi, . . . , mx}, with J2j — ^ (^^^ pl^t legends for a pictorial 
representation of the molecular groupings). Panel (a) depicts the instance M = N = 6, K = 3; the molecules are made of, from bottom to top, 
{1, 1, 4} modes (solid red line), {1, 2, 3} modes (dashed green line), and {2, 2, 2} modes (dotted blue line), respectively. Panel (b) illustrates 
the instance M = = 8, = 4; the molecules are made of, from bottom to top, {1,1,1,5} modes (solid red line), {1, 1, 2, 4} modes 
(dashed magenta line), {1,1,3,3} modes (dot-dashed green line), {1,2,2,3} modes (dot-dot-dashed cyan line), and {2,2,2,2} modes (dotted 
blue line), respectively. Panel (c) depicts the instance M = N = 15, = 5; the molecules are made of, from bottom to top, {1,1,1,1,11} 
modes (solid red line), {1, 1, 1, 3, 7} modes (dashed magenta line), {1, 2, 3, 4, 5} modes (dot-dashed green line), {2, 2, 2, 3, 6} modes (dot- 
dot-dashed cyan line), and {3, 3, 3, 3, 3} modes (dotted blue line), respectively (only 5 out of 30 possible molecular partitions are shown). 
Panel (d) illustrates a mixed-state instance: the total number of modes in the pure state is = 28 but only M = 24 of them are partitioned 
into K = 6 molecules (i.e. we are tracing over the extra 4 modes to be left with the multipartite molecular entanglement in a 24-mode 
fully symmetric mixed state); the molecules are made of, from bottom to top, {1,1,1,1,1, 19} modes (solid red line), {1,1,2,3,8,9} modes 
(dashed magenta line), {2, 2, 2, 6, 6, 6} modes (dot-dashed green line), {3, 3, 3, 3, 4, 8} modes (dot-dot-dashed cyan line), and {4, 4, 4, 4, 4, 4} 
modes (dotted blue line), respectively (only 5 out of 199 possible molecular partitions are shown). Notice how, for a fixed total number of 
modes and parties (molecules), the distributed multipartite entanglement increases with the ranking of the sorted list of partitions, i.e. it is 
bigger whenever the individual molecule sizes are more similar to each other. Notice also how the genuine multipartite entanglement in any 
partition increases unboundedly with the squeezing r ^ for pure states [Panels (a)-(c)], while it eventually saturates for mixed states [Panel 
(d)]. All the quantities plotted are dimensionless. 
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The standard form of multi-GLEMS corresponds to a K- reduced states of the molecules Mi and Mj (made of rrii 

and rrij modes respectively) of the original fully symmetric 



mode CM o"^^] , where the reduced CM of any two 

K:{mi,...,mK} . f 

modes i and j denotes a GLEMS in standard form 11571] . pa- 
rameterized by the local purity 1/a of mode i, the local purity 
1/6 of mode j, and the global purity 1/c of the two-mode re- 
duced state: 



/a ^+0 
a g- 
b 
\0 g- b 



(17) 



w ith g± = Wija- bf - (c - 1)^] [(g - bf - {cTm ^ 
^[{a + bf - (c - 1)2] [(a + bf - (c + 1)2] } / [4v^ . The 
three purities are in turn related to the determinants of the 



(AT) 



State cr c 

tion to construct the standard form CM o"^^] . 

K :|mi ,...,771x1 



via Eq. ([T2l) . A Mathematica instruc- 

of 

multi-GLEMS — obtained from unitary localization applied 
to fully symmetric states ^''^j^^ — given N, K, and 

a partition {mi, . . . , m^}, is given by Eq. (A. 10) in the Ap- 
pendix. To provide an illustrative example, the following CM 
describes a mixed 5-mode multi-GLEMS obtained by starting 
with a pure fully symmetric 20-mode Gaussian state, Eq. (fTOl) . 
with squeezing r = 1.0; discarding 5 modes; partitioning the 
remaining 15 modes into molecules, all of different size; and, 
finally, performing the unitary localization: 



~/(20) 

^5:{1,2,3,4,5} 



/ 1.871 1.658 1.834 1.956 2.044 \ 

1.871 -0.1587 -0.2152 -0.2691 -0.3217 

1.658 2.395 2.232 2.380 2.488 

-0.1587 2.395 -0.3535 -0.4421 -0.5287 

1.834 2.232 2.776 2.633 2.752 

-0.2152 -0.3535 2.776 -0.5996 -0.7169 

1.956 2.380 2.633 3.069 2.934 

-0.2691 -0.4421 -0.5996 3.069 -0.8966 

2.044 2.488 2.752 2.934 3.296 

\ -0.3217 -0.5287 -0.7169 -0.8966 3.296 / 



This state is completely non- symmetric (the reduced single- 
mode purities are all different) and it is characterized by a 
global purity /i ^ 0.30 and a genuine multipartite entangle- 
ment shared by the 5 modes equal to G^^^ ^ 0.014, as analyt- 
ically computable from Eq. ([T5]) [coded as Eq. (A.9)], given 
the properties of the original, un-localized state cr^^.^\ 2 3 4 5}- 
Thanks to the mechanism of unitary localization, we will 
now illustrate how the results of the previous Section yield 
in fact strong evidence that non- symmetric (mixed) Gaus- 
sian states in the form of multi-GLEMS obey the strong 
monogamy decomposition Eq. (|5]) and that G^^^ is a proper 
quantifier of genuine multipartite entanglement even for 
Gaussian states that do not satisfy the constraint of permu- 
tation invariance. To substantiate this statements, we review 
and put in perspective some of the results obtained in the pre- 
vious Section. The genuine i^T-partite entanglement shared by 
K modes of a multi-GLEMS increases with the mixedness 
of each individual mode (these quantities are monotonically 
increasing functions of the squeezing r), similarly to what 
happens for pure- state bipartite entanglement in general quan- 
tum systems; it decreases with the global mixedness of the 
state; it increases with increasing symmetry among the differ- 
ent modes (i.e. with increasing similarity between each local 
mixedness). Concerning the scale invariance of the genuine 
multipartite entanglement, it can be best understood by not- 
ing that the fully symmetric A^-mode state ^''^j^^ | and 



there-scaled (5A')-modeone ^^'^jsmi stuk} correspond to 
the same effective i^T-mode unitarily localized multi-GLEMS 
^'io{mi uik}' ®' ^P^^^ from the additional, irrelevant, 
single-mode thermal states which are obviously s-degenerate 
in the latter case. 

We conclude our discussion by reporting the results of 
extended numerical tests on strong monogamy, i.e. on 
the non-negativity of the genuine multipartite entanglement 
[Ineq. (|7])], as thoroughly illustrated and discussed in Fig. [71 
The obtained results can be equivalently regarded in terms 
of molecular block entanglement in fully symmetric (pure 
or mixed) Gaussian states, or entanglement shared by sin- 
gle modes in non-symmetric (pure or mixed) multi-GLEMS. 
Adopting the first picture for simplicity, we generated states 
of the form ^^^j^^ mx}' ^^^^^ parameters were ran- 
domized, as explained in detail in the caption of Fig. [71 The 
fact that in all instances we found a well defined, non-negative 
G^^^, provides good evidence in support of the conjecture 
that Eq. ilSh is a proper quantifier of genuine multipartite en- 
tanglement, and that entanglement distributes according to a 
strong monogamy law of the form Eq. ©, independently of 
symmetry constraints under permutation of the modes. A full 
analytic proof of this claim is however, to date, available only 
in the special case of the multipartite entanglement shared by 
single modes in states cr^^l-^ ^| f2§\, and it is still lacking 
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FIG. 7: (Color online) Plot, as a function of the squeezing r, of 
the genuine multipartite entanglement G^^^ in 10000 random fully 
symmetric A^-mode Gaussian states. Each point corresponds to the 
K-partite entanglement among K molecules Mi, . . . , Mk, each 
comprising mj modes, with J2f^i — ^^^^ instance, the 

real parameter r is randomized between and 2, the integer N (to- 
tal number of modes of a pure fully symmetric Gaussian state with 
squeezing r) is randomized between 2 and 100, the integer M (ef- 
fective number of modes partitioned in molecules) is randomized be- 
tween 2 and N, and the integer K (number of molecules) is ran- 
domized between 2 and M. For any M < the multipartite en- 
tanglement is being computed in a mixed state (as (N — M) modes 
are traced out); furthermore, for any < M the molecules have 
in general different sizes. The partition {mi, . . . , rriK} is chosen at 
random as well among all the possible partitions. Upon application 
of the unitary localization, each point equivalently corresponds to the 
genuine i^-partite entanglement of generally mixed, generally non- 
symmetric random K-mode Gaussian states belonging to the class of 
multi-GLEMS. In up to two million states (not shown here) randomly 
generated according to this prescription, with up to 10^, all values 
of Gr^^ have been found non-negative. This numerical evidence sup- 
ports the conjecture that the strong monogamy decomposition holds 
for all Gaussian states beyond symmetry constraints under permuta- 
tion of the modes, and that the positive-defined residual contangle is 
a proper quantifier of genuine multipartite entanglement for general 
states on harmonic lattices. The solid line in the figure corresponds 
to the instance M = = 2, = 2, mi — m2 — 1, i.e. to 
the bipartite entanglement in two-mode squeezed states. As the gen- 
uine multipartite entanglement decreases with A^, such a bipartite 
entanglement stands as an absolute upper limit for any (A^ > 2)- 
mode instance. This bound can be only attained for A^ = M = 2s, 
K = 2, mi = m2 = s (with s an integer scale factor) by virtue of 
the scale invariance discussed in Sec. IV CI All the quantities plotted 
are dimensionless. 



in the general case of distributed molecular entanglement. 



VII. CONCLUDING REMARKS 

In the present work we have discussed the characteriza- 
tion and quantification of genuine multipartite entanglement 
in fully inseparable Gaussian states of multimode, bosonic CV 
systems. We adopted the framework introduced in Ref. ll28ll 
in order to establish some generalizations of the conventional 
monogamy inequalities on distributed entanglement and to de- 



fine a proper and computable quantifier of genuine multipar- 
tite entanglement. Our main aim has been to verify whether 
general Gaussian states, not endowed in principle with special 
symmetries with respect to mode permutation and/or not nec- 
essarily partitioned into parties each owning a single mode, 
obey strong monogamy constraints, i.e. that hold simultane- 
ously on bipartite and genuine multipartite entanglement. This 
type of analysis is motivated by conceptual questions as well 
as by the practical utility of multipartite Gaussian resources 
for CV quantum information and communication. Such a 
study is thus of value in order to single out the most efficient 
ways to engineer and partition Gaussian states yielding opti- 
mal performances in quantum information tasks. 

In this paper we did not present a comprehensive, rigorous, 
and conclusive analysis on the structure of multipartite entan- 
glement sharing and on the quantification of genuine multi- 
partite entanglement in general Gaussian states. On the other 
hand, we obtained a series of partial results that, in our opin- 
ion, points to a unified framework and realizes concrete steps 
towards a positive answer to the questions posed. To summa- 
rize, we considered a class of Gaussian states a^J^] . , 

' K:{mi,...,mK } 

with the following properties. They are M-mode fully sym- 
metric mixed Gaussian states (i.e. invariant under permutation 
of any two modes) obtained from N-mode fully symmetric 
pure Gaussian states by discarding N — M modes. These 
states are partitioned into K parties, or molecules, respec- 
tively formed by mi, ... , tuk modes (with Y^f^i = M. 
We derived an analytical formula — and we included a Math- 
ematica code to evaluate it numerically in arbitrary cases — 
for the genuine multipartite entanglement (quantified by the 
residual contangle G^^^ emerging from the strong monogamy 
decomposition) shared by the K parties. We characterized 
various properties of this quantity that we list here synopti- 
cally in order to provide a final comprehensive summary: 

(i) The genuine multipartite entanglement among K 
"molecules" of modes with fixed sizes ruj and assigned 
partition decreases if the global number of modes N of 
the original pure state is increased (resulting in a more 
mixed state of the M modes on which the partition is 
imposed). 

(ii) The genuine multipartite entanglement is invariant un- 
der global re-scaling of the system, i.e. varying by the 
same scale factor both the number of modes N and the 
individual sizes mj of each molecule, at fixed number 
K of molecules. 

(iii) For a fixed N and fixed number M of non-discarded 
modes, the genuine multipartite entanglement decreases 
with increasing number K of molecules in which the M 
modes are partitioned. 

(iv) For fixed N, M, and K, the genuine multipartite en- 
tanglement exhibits a hierarchical structure with respect 
to the individual sizes of the molecules rrij . It appears 
to increase with increasing ranking (order index) of the 
sorted list of all possible partitions of M modes into K 
molecules. Roughly speaking, it increases as the sizes 
of individual molecules get closer. 
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(v) The genuine multipartite entanglement always in- 
creases with the average squeezing r for fixed values of 
the other parameters. For any multi-partition of glob- 
ally pure states (M = N) the multipartite entanglement 
diverges in the limit of infinite squeezing, while for any 
multi-partition of globally mixed states (M < N) it 
saturates to a finite value. In all states of the considered 
class there is thus a promiscuous entanglement distribu- 
tion, i.e. a coexistence between bipartite and any form 
of multipartite entanglement, with all contributions be- 
ing increasing functions of each other. 

(vi) The genuine multipartite entanglement among K 
molecules with, generally different, sizes mi , . . . , mx 
in pure or mixed fully symmetric Gaussian states, is 
unitarily localizable, i.e. it is equal to the genuine multi- 
partite entanglement among K modes in pure or mixed, 
generally non- symmetric, Gaussian states with a well 
defined symplectic spectrum characterized by all sym- 
plectic eigenvalues equal to unity but one. 

(vii) The genuine multipartite entanglement is found to be 
non-negative in all the instances analyzed, and af- 
ter extensive numerical tests on ample classes of ran- 
domly generated permutation-invariant CMs. This re- 
sult provides evidence that the strong monogamy con- 
straint holds true both for the bipartite and the genuine 
multipartite entanglement distributed among blocks of 
modes of arbitrary size in fully symmetric Gaussian 
states, and among single modes in the subset of non- 
symmetric Gaussian states constituted by mixed states 
of partial minimum uncertainty. 

Property (vi) stems from a more general result, of indepen- 
dent interest, demonstrated in this paper, namely the unitary 
localization of bipartite and multipartite correlations in multi- 
symmetric Gaussian states (i.e. states that are locally invariant 
under permutation of modes within each molecule, but are in 
general not invariant under permutation of two modes taken 
from different molecules), which extends an analogous result 
proved previously in the bipartite instance |32]. 

The above findings, together with the operational interpre- 
tation established in Ref. |28] that relates G^^^ to the optimal 
fidelity of CV teleportation networks 1 10, 23] employing fully 
symmetric Gaussian resources, yield further support to the use 
of the generalized residual contangle G^^^ for the quantifica- 
tion of genuine multipartite entanglement in Gaussian states. 

Whether the strong monogamy of bipartite and genuine 
multipartite entanglement is a property that applies in gen- 
eral to all quantum states remains an open question. From the 
present results in the Gaussian CV scenario, and in particular 
from point (vii) above, one is lead to conclude that the con- 
straints on entanglement imposed by the strong monogamy 
have a wider domain of validity than de Finetti-type bounds 
1 47, 59] which apply only to fully permutation-invariant 
states, even though the monogamy-based bounds carry the 
drawback of being crucially dependent on the actual mea- 
sure of entanglement that is chosen. An investigation of the 
strong monogamy decomposition and an evaluation of the cor- 
responding quantifier of genuine multipartite entanglement on 



states of finite-dimensional systems would provide a relevant 
test ground in order to verify the general validity of the prop- 
erty of strong monogamy of distributed quantum correlations. 
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APPENDIX: MATHEMATICA CODE 

This appendix lists some of the routines we defined in Wol- 
fram Mathematica 6.0 and used for the evaluation of the 
genuine multipartite entanglement shared by modes and/or 
molecules in fully symmetric Gaussian states. A brief com- 
ment accompanies each instruction; more details can be re- 
trieved in the main text, in the specific part in which the cor- 
responding instruction is referred to. The following two pack- 
ages need to be loaded in the Mathematica kernel: 

<< Linear AlgebraMattixManipulatiori 
<< Combinatorica 



(A.l) Definition of contangle in terms of the squared deter- 
minant (P involved in the computation of the Gaussian 
entanglement measures [see Eq. ©] : 

contangle[J2_] := (ArcSinh[Sqrt[-l -h ^/2]])2; 



(A. 2) Covariance matrix of fully symmetric pure Gaussian 
states of n modes: 

FulSymCM[n_] := BlockMatrix[Table[ 

(V«n-l))(<*'r'<»;!^,)'- 

{kl^2-\- 9{-2 + n)-n,k2^ 

Sqrt[(-l + 62) (62^2 -(n- 2)2)]})], 
{i,l,n},{j,l,n}]]; 

(A. 3) Relation between single-mode determinant iP' and 
squeezing r of pure, fully symmetric n-mode Gaussian 
states, as a function of n: 

mybsq[n_] := 4^(2 -|- (-2 -h n)n -h 2(-l -h n)Cosh[4r]); 

(A.4) Reduced determinant of a block of / modes out of an 
n-mode fully symmetric pure Gaussian state: 
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ghredet[n_, L] := 

(l/(n - 1)) * {l{n - l)h^ -(I- l)(n - Z - 1)) ; 
ghd[n_, L] := 

Re&ie[Simplify[gliredet[n, l]Lh Sqrt[mybsq[n]]], 
r > 0]//Simplify; 



(A. 5) Gaussian entanglement measure (squared determinant 
d^) for a two-mode GLEMS with local purities 1/a and 
1/b and global purity 1 / c: 

d2glems[a_, c_] := 

If [c > V-l + a^ + ft^, l,If[c < Sqrt[ 

(2 (a^ + b') + (a^ - 6^) V Abs [a^ - b^] 

V((a^-6^f + 8(a^ + 6^)))],g^, 

gi^(-a4-64- V((-l + «-&-c)(l + a-6-c) 
(-1 + a + 6 - c)(l + a + 6 - c))v'((-l + a - 6 + c) 
(1 + a - 6 + c)(-l + a + 6 + c)(l + a + 6 + c))- 

(-1 + c^)' + 26^ (1 + c^) + 2a2 (1 + 6^ + c^))]] ; 

(A. 6) Gaussian entanglement measure (squared determinant 
d'^) for a (z|j)-mode bipartition in a n-mode pure fully 
symmetric Gaussian state: 

d2ijn[n_,i_,i_] := 

Refine[(((Refine[(d2glems[Sqrt[ghredet[n, i]] , 
Sqrt[ghredet[n, jj] , Sqrt[ 
ghredet[n, z + j]]]//Simplify), b > 1])/. 
6 ^ Sqrt[mybsq[n]])//Simplify), r > 0]; 

(A.7) Genuine /c-partite entanglement shared by k individ- 
ual modes in a n-mode pure fully symmetric state, as 
a function of /c, n, and the squeezing r (as computed in 

ii): 

GEnt[A:_, n_, r_] := 

Sum [(-l)^ ArcSinh [ {2^/-l - j -\- kSinh[2r]) / 
Binomial[-l -h {i, 0, - 2}]; 

(A. 8) Sorted list of all the possible partitions of n modes into 
k molecules: 



KSoitedPartitions[n_, := Modulef 
{tabpart^ tabpartsort^ u}^ tabpart = Partilions[n]; 
tabpartsort = {};Do\H\]^ngth[tabpart[[u]]] == fc, 
AppendTo[tabpartsort^ {{Sort[tabpart[[u]]]}}]]^ 
{n, lj&ngfii[tabpart] }] ; 
Retum[MatrixFoim[Sort[^aA>parrsoit]]]] ; 



(A.9) Genuine multipartite entanglement shared by k 
molecules made of mpart = {mi, rrik} modes 

respectively mj < , in a global n-mode pure 

fully symmetric Gaussian state (recursive definition): 
GMolEnt[n_, mparr_] := 
contangle[d2ijn[n, mpar/[[l]], 
Siim[mparr[[i]], {i, 2, k}]]] — Sum[ 
Sum[GMolEnt[n, m, KSubsets [mp/2ft, m] [[Z]]], 
{Z,Binomid[fc — l,m — l]}],{m, fc— 1,2,-1}]; 



(A. 10) Covariance matrix of the equivalent /c-mode state 
(multi-GLEMS) associated, via unitary localization, 
to a multi-partition of n modes into k molecules 
made of mpart = {mi, . . ., rrik} modes respectively 

(Ej rrij < nj : 

UiiiLocCM[n_, mpart_] := 
Refine [(Refine [(BlockMatrix [(Table [If [i == j, ai, 
If < ^ij^^oA] ^ k}, (j, k}])] /. 
Join [Table [ai Sqrt[ghd[n, mpa/t[[i]]]] 
IdentityMatrix[2], {i, fc}], 
Flatten [Table [cij DiagonaIMattix[ 

{4;^ W ((- (-1 + atpj) 2 + K - aj) 2) 
(-(l + aip,)2 + (ai-a,)2)) + 

(- (1 + + ^)) - 

^{{-{-l + ai^y + ioi + aj)^) 
(-(l+a,p,) 2 + + 2)))}/. 
{tti ^ Sqrt[ghd[n,mparr[[i]]]],aj 
Sqn[ghd[n,mpart[\j]]]],aipj Sqrt[ 
ghd[n,mpaft[[i]] +mpaft[|>]]]]}], 
{i,fc},{i,fc}]]]),r>0]// 
Simplify), r > 0]//Simplify; 
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